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Abstract: This paper presents a new method for solving system oj higher -order linear differential equations (HLDEs) with constant coefficients; 
however the same idea can be extended to variable coefficients. Using the basic concept of inverse of matrix and variation of parameters, we develop a 
new method to solve system of HLDEs. Proposed method works for any right hand side function, so-called vector forcing function f (x) of given system of 
HLDEs. Selected examples are presented using proposed method to show the efficiency. 
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I INTRODUCTION 

The systems of higher-order linear differential equations (HLDEs) have been vigorously pursued by many researchers and engineers 
and developed many different method to solve the system, for example, see [1-10]. Naturally, the systems of HLDEs arise in many 
applications of nuclear reactors, multi-body systems, vibrating wires in magnetic fields, models of electrical circuits, mechanical 
systems, robotic modeling, and diffusion processes etc. 

We consider the following type of systems of n linear differential equations of order m > 0, 



d m cl 
A m T^^W^" * + A —u(x) + \u(x) = f(x\ 
ax ax 



(1) 



where, for i — 0, ... , n, A i are coefficient matrices of order nxn , f(x) = (f^x), ... , f n (x)) T and u(x) = (u x (x), ...,U n (x)) T are 
an n- dimensional vector forcing function and unknown vector respectively. If the lading coefficient A m is non-singular, then the system 
of HLDEs (1) is called as of first kind, otherwise it is called second kind. In this paper, we present a new method to solve the system (1) 
of first kind with constant coefficients; however the same idea can be extended to variable coefficients. For obtaining the general 
solution of given system, we use the basic concepts of inverse of matrix and variation of parameters formula. 

Rest of the paper is organized as follows: In Section II, we present the proposed method to solve system of HLDEs and selected 
examples are presented in Section III to show the efficiency of proposed method. 



II A NEW METHOD FOR SYSTEM OF HIGHER-ORDER LINEAR DIFFERENTIAL 

EQUATIONS 

Recall the system of HLDEs given in Section I 
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d m d 
L(u(x)) = A m -— u(x) + • • • + \ —u(x) + \u(x) = /(*), 
ax ax 



(2) 



where L — A m 1 h A x h Aq, is nxn matrix differential operator. Now the following theorem presents the algorithm to 

dx m dx 

compute the general solution of system (2), with given set of fundamental system of determinant of matrix differential operator I. 
Theorem 1 : Given a fundamental system {^(x), . . . , V mn (x)} , of T = det(L) , then the system of higher -order LDEs 



has following solution 



d m d 
A m ~T~^~ W ( X ) + • • • + A — u(x) + AqU (x) = f(x), 
dx dx 



(3) 



u(x) - 



det(Z})f> y (*)J 



det(L*)f> 7 (x)j" 



detg^(x) 
det(W) 

det(W,)/;.(x) 
det(W) 



</.v 



(4) 



where Wis the Wronskian matrix of {^(x), . . . , V mn (x)} and W i obtained from Why replacing the i-th column by mn-th unit vector; 
and det(Z^) denotes the determinant of L after removing i-row and k-th column. 

Proof: Since the leading coefficient A m of system (3) is non-singular matrix, the inverse of L exists with det(I) = T as scalar differential 
operator of order mn. If {^(x), ... , V mn (x) } is fundamental system of T. Now the solution of system (3) can be obtained as 

= Adj(L)y(x) , where y(x) = (y^x), ... , y n (x)) T is solution obtained from Ty t (x) = /. (x) computed using the classical 

formulation variation of parameters as follows. The scalar differential equation Ty(^x} = /(-^) can be reformulated as system of first 

order linear differential equation, say y\x) = My(x) + f(x) , where M is companion matrix. Now the solution of first order system is 

obtained as y(x) = wj W~ l f(x) dx and the solution Ty (x) = / (x) is the first row of y(x) . 

The solution (4) satisfies given system of HLDEs, as follows: 

det(W.)/;.(x) 



L(u(x)) — L 



det(Z})|>,(jc)J 



det(L»)f>.(x)j" 



det(W) 



dx 



fet<Wj)fiW 



7=1 



det(W) 

^(-l) 1+1 det(L;> .» (-l) 1+ "det(Zf) N 
det(4 ) • • • (-l) n+n det(L* ) 



dx 



> v.(x)\ dx 

pi J J det(W) 



det(W. )/„(*) 
det(W) 



=(/„ ...,/„)"=/ 



Therefore, the solution (4) is general solution system (3). 

In following section, we present selected numerical examples (system of HLDEs with constant coefficients and variable coefficients) 
using proposed method presented in Theorem 1 . 
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III NUMERICAL EXAMPLES 



In this section, we present couple of examples for system of linear differential equations with constant coefficients (Example 1) and 
variable coefficients (Example 2) respectively, to show the efficiency of proposed method in Theorem 1 . 

Example 1 : Consider the following system of linear differential equations of order two with constant coefficients 

d 2 



u x (x) + u x (x) - 2u 2 (x) = sin(x), 
dx 

d 2 , x d . . _ 3x 

— yU 2 (x) u 2 (x) = e . 

dx dx 



(5) 



The system (5) can be written in matrix notations as L(u(x)) = f(x), with L = — - + \ h Aq , where 



dx dx 



4 



0 1 



'0 0^ 
0 -1 



> A> 



0 0 



and u(x) - 



r u x {x)^ 

\U 2 (x)j 



.fix): 



sin(x) 

-3x 



V 



If we denote D = , then L ~- 

dx 



r D 2 +l -2 " 



Following procedure in Theorem 1 , we have 



0 D 2 -D y 
Ty(x) = D 4 y(x)-D 3 y(x) + D 2 y(x)-Dy(x) = f(x) , 
here y(x) = (y l (x), y 2 (x)) T . Now from Theorem 1 , we have the general solution of given system (5) as follows 

1 



u(x) ■ 



(~u x (x) 
K u 2 (x)j 



c, cos(x) + c 9 sin(x) — cos(x) + 2c, + c A e x H e x 

1 2 2 3 60 

1 -3x 

— e 
12 



where c 1 ,c 2 ,c 3 and c 4 are arbitrary constants. 

Example 2: Consider the following system of linear differential equations of order one with variable coefficients 



— u l (x) — x—u 2 (x) + u l (x) — (l + x)u 2 =0, 
dx dx 

d s \ /N 2x 

— u 2 (x) + u 2 (x) = e . 
dx 



(6) 



The matrix representation of system (5) is L(u(x)) = \ u(x) + AqU(x) — f(x) , where 

dx 



A = 



1 -x 
v 0 \j 



1 -l-x 
0 1 



and u(x) ■ 



r u x (x)^ 

yU 2 {x)j 



fix)-- 



Matrix differential operator L = 
given system (6) as follows 



D + l 



-xD-l-x^ 
D + l 



, where D = . Following Theorem 1 , we have the general solution of 

dx 
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u(x) = 



u y {x) 
K u 2 (x) j 



X 1 

c x e~ x + 2c 0 xe~ x + c 0 e~ x + -e 2x + -e 2x 
1 2 2 3 9 

c 1 e~ x + -e 2x 
2 3 



where c x and c 2 are arbitrary constants. 



IV CONCLUSION 



In this paper, we have presented a new method for solving system of higher order linear differential equations. This proposed method 
works for any vector forcing function of given system of HLDEs. The proposed method is developed using the basic concept of inverse 
of matrix and variation of parameters. Couple of examples (systems with constants coefficients and variable coefficients) are presented 
using proposed method to show the efficiency. 
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